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NOMENCLATURE 


B  ship  beam 

Cn  dimensionless  ice  strength 
D  propeller  diameter 
/  apparent  ice-hull  friction  factor 
/j  ice-hull  friction  factor 

F  Froude  number  based  on  ice  thickness:  V'/Ve// 

n  o  i 

g  acceleration  of  gravity 

h.  ice  thickness 
1 

h  snow  thickness 

S 

/i  apparent  advance  coefficient:  VjuD 
Ky  torque  coefficient 
Kt  thrust  coefficient 

average  rate  of  propeller  rotation 
N  normal  pressure  on  ice  sample  during  friction  test 
PD  delivered  power 

average  propeller  torque 
R  net  resistance  in  ice:  R,  -  R 
R  total  resistance  in  level  ice 

it 

R  resistance  in  ice-free  water 

OW 

R  predicted  total  resistance 
K,  nondimensional  resistance 
S  sampling  rate  of  data  acquisition  (milliseconds) 
t  thrust  deduction  coefficient 
T  ship  draft;  also  average  tangential  force  on  ice 
sample  during  friction  test 
Ta  average  propeller  thrust 

Tj  t  average  load  exerted  on  either  load  cell  during  friction  tests 
Tq  initial  tangential  force  due  to  ice  adhesion 
V ,  v  ship  speed 

IV  weight  applied  to  ice  sample  during  friction  tests 
IV  weight  of  ice  sample  and  sample  holder 
IV  total  weight:  !V+  W'^ 
v  Poisson's  ratio  of  ice  (~  0.3) 

X  geometric  model  scale 
p  water  density 
p(  ice  density 
y  water  specific  weight 
oi  ice  flexural  strength 
x  ice  friction  stress 

x  ice  friction  stress  ascribable  to  adhesion 

O 


V 


Model  Tests  in  Ice  of  a  Canadian  Coast  Guard 
R-class  Icebreaker — High  Friction  Model 

JEAN-CLAUDE  TATINCLAUX  AND  CARL  R.  MARTINSON 

INTRODUCTION 

At  the  1 6th  International  Towing  Tank  Conference  (ITTC),  Leningrad,  USSR,  1981,  the 
Committee  on  Performance  of  Ships  in  Ice-Covered  Waters  (Ice  Committee,  for  short)  in¬ 
itiated  an  international  cooperative  research  program.  Under  this  cooperative  program,  the 
participating  ice-testing  facilities  were  to  conduct  a  specified  test  program  with  the  same 
icebreaker  model  according  to  their  own  testing  and  data  analysis  methods  (ITTC  1981). 
Two  1:20  scale  models  and  two  1>0  scale  models  of  the  Canadian  Coast  Guard  R-class  ice¬ 
breaker  were  specially  constructed  by  the  National  Research  Council  of  Canada  and  cir¬ 
culated  among  the  participating  facilities.  The  smaller  1:40  scale  models  were  to  be  tested 
in  resistance  only  at 

•  Arctic  and  Antarctic  Research  Institute  (AARI),  Leningrad.  USSR. 

•  Norwegian  Technical  Institute  (NTI),  Trondheim,  Norway. 

while  the  larger  models  were  to  be  tested  in  resistance  and  propulsion  at 

•  CRREL,  Hanover,  New  Hampshire,  USA. 

•  Hamburgische  Schiffbau  Versuchsanstalt  (HSVA),  Hamburg,  FRG. 

•  Tapan  Ship  Research  Institute  (JSRI),  Tokyo,  Japan. 

•  Institute  for  Marine  Dynamics/National  Research  Council  (NRCC), 

St-Johns,  Newfoundland,  Canada. 

•  Wartsila  Arctic  Research  Center  (WARC),  Helsinki,  Finland. 

Resistance  and  propulsion  tests  in  ice-free  water  had  been  made  on  an  earlier  1 :40  scale 
model  and  previously  reported  by  the  National  Research  Council  of  Canada  (Murdev  1980). 

The  results  of  the  tests  on  the  original  ship  models  were  reported  at  the  1984  ITTC  in 
Goteborg,  Sweden  (ITTC  1984).  Comparison  of  the  test  results  from  the  various  facilities 
showed  that  under  nominally  identical  conditions  of  ice  thickness,  ice  strength  and  ship- 
model  speed,  the  ice  resistance  and  the  propeller  thrust  and  torque  of  the  1:20  scale  model 
were  within  25%  of  one  another.  On  the  other  hand,  when  these  test  results  were  ex¬ 
trapolated  to  full-scale  conditions,  the  predicted  resistance,  thrust  and  torque  were  signifi¬ 
cantly  below  available  full-scale  trial  measurements  (Edwards  et  al.  1981,  Michailidis  and 
Murdey  1981).  The  discrepancies  between  predicted  performance  and  full-scale  measure¬ 
ments  were  attributed  to  the  ice  friction  factor  of  the  model  hull,  measured  at  0.04  in  the 
average,  being  much  lower  than  the  estimated  value  of  0.1  for  a  new  icebreaker  hull.  The 
ITTC  Ice  Committee  decided  to  repeat  the  resistance  and  propulsion  tests  with  roughened 
1:20  scale  models.  WARC  agreed  to  treat  the  model  hulls  to  achieve  a  friction  factor  of 


approximately  0.1  and  to  prepare  a  friction  test  board  in  a  similar  fashion.  Nippon  Kokan 
Tsu  Laboratories  (NKK),  Tsu  City,  Japan,  who  had  recently  inaugurated  an  ice  towing  tank, 
joined  the  original  facilities  in  this  new  phase  of  the  cooperative  test  program.  The  results 
of  this  second  series  of  tests  were  reported  at  the  1 8th  1TTC,  October  1987,  Kobe,  Japan  (1TTC 
1987). 

The  results  of  the  test  series  conducted  at  CRREL  on  the  smooth  model  have  been 
reported  earlier  (Tatinclaux  1984).  This  report  presents  the  results  of  the  resistance  and 
propulsion  tests  performed  at  CRREL  with  the  roughened  model  and  compares  the  test 
results  obtained  with  the  two  models  at  the  facilities  involved  in  the  cooperative  test 
program. 

SHIP  CHARACTERISTICS  AND  TEST  CONDITIONS 
Ship  characteristics 

The  Canadian  R-class  icebreaker  was  designed  to  operate  continuously  in  1-m-thick  level 
ice.  The  ship  has  a  displacement  of  approximately  8,000  tons  at  a  midship  draft  of  6.9  m  and 
is  propelled  by  twin  fixed-pitch  propellers  with  a  total  shaft  power  of  1 1,000  kW.  The  main 
hull  and  propeller  characteristics  at  full  and  model  scales  are  listed  in  Table  1.  A  photo¬ 
graphs  of  the  model  is  shown  in  Figure  1.  Three  R-class  icebreakers  are  in  operation,  the 
CCGS  Rinlbson,  theCCGS  Franklin  and  the  CCGS  Dc>  Groscillcr s.  Full-scale  trials  have  been 
conducted  and  reported  by  Edwards  et.  al.  (1981)  for  the  Radbaon  and  bv  Michailidis  and 
Murdev  (1981)  for  the  Franklin. 


Table  1.  Mean  characteristics  of  R-class  icebreaker. 
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Test  conditions 

The  set  of  conditions  for  the  resistance  and  propulsion  tests  to  be  performed  with  the  1 :20 
scale  model,  as  agreed  to  by  the  ITTC  Ice  Committee  members,  is  given  in  Table  2. 


Table  2.  Test  conditions. 

Full  settle  Mo, Icl  settle 

Ice  thickness,  h  45  and  70  cm  22.5  and  35  mm 

Ice  strength,  <r  400  and  800  kPa  20and40kPa 

Ship  speed,  V  0.5to5.0m/s  0.1 1  to  1.12  m/s 

Range  ot  Fn  0.2  to  2.4 

Range  of  C  n()  to  ISO 


ICE  FRICTION  FACTOR 


As  mentioned  previously,  WARC  had  treated  both  the  model  hull  and  a  special  friction 
plate  to  achieve  an  ice  friction  factor  of  approximately  0.1.  Friction  tests  were  conducted 
using  the  plate  mounted  on  a  recently  built  friction  table.  An  ice  sample,  1 3.5  bv  1 3.5  cm  in 
plan  dimensions,  was  inserted  into  the  fixed  sample  holder  that  was  connected  to  two  load 
cells.  A  weight,  W,  selected  to  exert  the  required  normal  pressure,  N,  was  placed  on  the 
sample  holder.  Both  load  cells  were  pre-tensioned  and  the  friction  table  was  set  in  motion 


at  escribed  speed.  In  most  of  the  friction  tests,  the  table  traveled  back  and  forth  twice  to 
account  for  possible  minor  misalignment  of  the  friction  table  with  respect  to  the  horizontal 
plane.  The  average  frictional  force,  T,  was  measured  as  the  mean  of  the  average  forces 
recorded  by  the  two  load  cells,  T  =  (T;  +  T2)/2,  over  a  full  cycle.  The  apparent  friction 
coefficient,  / ,  was  calculated  as  the  ratio 


/.  =  T/(W+W ) 

'a  S 


(1) 


where  W  was  the  weight  of  the  ice  sample  and  sample  holder.  The  friction  test  conditions 
and  results  are  listed  in  Table  3.  All  tests  were  made  at  the  travel  speed  of  approximately 


Table  3.  Results  of  ice  friction  tests. 

Apparent  friction 


Normal  load 
W,  (N) 

Normal  pressure 
N  (1 (Pa) 

Tang,  force 
TIN) 

factor 

iJxlOO) 

69 

3.8 

4.9 

7.11 

69 

3.8 

5.2 

7.51 

69 

3.8 

5.7 

8.20 

69 

3.8 

5.7 

8.28 

69 

3.8 

4.6 

6.57 

69 

3.8 

4.7 

6.83 

69 

3.8 

5.0 

7.24 

69 

3.8 

5.2 

7.44 

36.5 

2.0 

2.7 

7.31 

36.5 

2.0 

2.9 

7.96 

36.5 

2.0 

3.1 

8.40 

36.5 

2.0 

3.2 

8.64 

16.4 

0.9 

1.4 

8.72 

16.4 

0.9 

1.6 

9.81 

16.4 

0.9 

1.7 

10.54 

16.4 

0.9 

1.8 

10.76 

16.4 

0.9 

1.7 

10.16 

16.4 

0.9 

1.8 

11.01 

16.4 

0.9 

1.7 

10.37 

16.4 

0.9 

1.9 

11.34 

Averages 

69.3 

3.8 

5.1 

7.4±0.8 

36.5 

2.0 

3.0 

8.1±1.3 

16.4 

0.9 

1.7 

10.413.4 

10  cm/s  and  under  wet  conditions  by 
pouring  a  thin  film  of  water  over  the 
test  plate.  No  friction  tests  were  made 
directly  on  the  model  hull.  The  aver¬ 
age  value  of  /  for  all  the  tests  per¬ 
formed  was  found  to  be  0.086,  with  a 
standard  deviation  of  0.026. 

The  friction  test  results  are  shown 
in  Figure  2a  as  /  versus  normal  pres¬ 
sure,  N.  It  can  be  seen  that  the  appar¬ 
ent  friction  factor  is  decreasing  with 
increasing  normal  pressure  or  load. 
Similar  behavior  of  f  with  N  had  been 
observed  in  a  previous  study  by  For- 
land  and  Tatinclaux  (1985)  who  at¬ 
tributed  this  phenomenon  to  the  exis¬ 
tence  of  an  adhesion  force  between 
the  ice  sample  and  the  test  surface,  the 
origin  of  which  is  not  yet  fully  under¬ 
stood.  This  adhesion  force  was  con¬ 
sidered  to  give  rise  to  an  additional 
tangential  stress,  T  ,  at  the  ice/ plate 
intr  rface,  the  magnitude  of  which  had 
been  found  to  increase  with  decreas¬ 
ing  ice  toughness  or  icestrength.  Thus, 
the  total  frictional  stress  can  be  ex¬ 
pressed  as 


T=/jN  +  To  (2) 

where  /.  is  the  actual  friction  factor.  Indeed,  when  t  is  plotted  against  N  as  in  Figure  2  it  is 
clear  that  X  is  a  linear  function  of  N.  Linear  regression  analysis  of  the  test  data  yielded  to  = 
66  Pa  and  f  =  0.071  compared  to  the  average  apparent  friction  factor  of  0.086,  that  is 

t(kPa)  =  0.066  +  0.071  N(kPa) .  (2a) 
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a.  Apparent  friction  factor  vs  normal  pro-  />.  Tangential  stress  vs  normal  pressure, 

sure. 

Figure!.  Results  of  friction  tests. 

RESISTANCE  TESTS  WITH  ROUGHENED  MODEL 

Experimental  conditions  and  test  results 

In  the  resistance  tests,  the  ship  model  was  connected  to  the  towing  carriage  by  a  rigid  7.6- 
ern  (3-in.)  diameter  towing  post  that  could  slide  vertically  in  a  linear  ball-bearing.  The  tow 
post  was  attached  to  a  force  block,  mounted  on  a  double-gimbal  that  was  bolted  to  the 
bottom  of  the  model.  The  ship  model  was  thus  free  to  heave,  pitch  and  roll  but  was  totally 
restrained  in  surge  and  sway.  It  was  restricted  n  yaw  by  a  fork  that  was  attached  to  the 
carriage  and  straddled  a  vertical  rod  at  the  stern  of  the  model.  The  ship  model  in  the  trim 
tank  of  the  ice  towing  tank  is  shown  in  Figure  I. 

For  each  ice  sheet,  three  tests  were  made  at  different  velocities,  each  over  a  distance  of 
about  1 0  m  or  two  model  lengths.  After  the  first  and  second  tests,  the  model  was  backed  into 
the  previously  broken  channel  a  distance  sufficient  toallowit  to  reach  steady  velocity  before 
it  entered  the  unbroken  ice.  Before  each  test  series,  the  ice  thickness  and  flexural  strength 
were  measured  at  several  locations  along  the  tank,  the  latter  by  the  in-situ  cantilever  beam 
method.  The  modulus,  or  characteristic  length  of  the  ice  sheet,  was  not  measured.  However, 
for  the  range  of  ice  strengths  used  in  the  study,  past  experience  has  shown  that  the 
characteristic  length  of  the  ice  sheet  is  approximately  10  times  its  thickness. 

The  results  of  the  resistance  tests  are  listed  in  both  dimensional  and  nondimensional 
form  in  Table  4.  The  model  resistance  in  clear  water,  R  ,  was  calculated  as 

OW 

R  =  14.53  V]MT  (3) 

where  R  is  expressed  in  Newtons  and  V  is  the  model  speed  in  meters  per  second .  Equation 
3  was  obtained  by  regression  analysis  of  the  resistance  test  results  bv  Murdey  (1980)  with 
a  correlation  coefficient  r  =  0.997.  Tests  at  JSRI  (ITTC,  1 987)  have  shown  negligible  effect  of 
the  hull  roughness  on  R  .  The  net  ice  resistance  was  then  calculated  as  R  =  R  -R  .  The 
clear  water  resistance  is  usually  very  small  as  compared  to  the  ice  resistance  and  well  with¬ 
in  the  range  of  uncertainty  of  R  .  Only  at  relatively  high  speed  in  thin  ice  does  R  become 
a  significant  component  of  R  . 

Data  analysis 

Because  of  the  lack  of  a  satisfactory  analytical  expression  for  ship  resistance  in  ice,  it  is 
customary  to  fit  empirical  relationships  to  the  nondimensional  data  using  regression 
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Table  4.  Results  of  resistance  tests  in  level  ice. 


Dimensional  data  Nondimensional  data 


Test  no. 

h 

/ 

(cm) 

a, 

(k  Pa) 

V 

(ml s) 

R„ 

<N) 

R 

!>tl» 

(N) 

R, 

(N) 

C„ 

F 

n 

VK 

R/yBh; 

101 

22 

28 

0.11 

23 

0.2 

22.8 

130 

0.24 

4.99 

4.95 

102 

17 

30 

0.34 

27 

1.7 

25.3 

180 

0.83 

9.82 

9.19 

103 

20 

20 

0.57 

40 

4.8 

35.2 

102 

1.29 

10.51 

9.25 

111 

24 

24 

0.78 

67 

8.9 

58.1 

102 

1.61 

12.22 

10.60 

112 

18 

25 

1.02 

67 

15.1 

51.9 

142 

2.43 

21.73 

16.83 

113 

23 

26 

1.23 

97 

21.8 

75.2 

115 

2.59 

19.27 

14.93 

121 

26 

21 

0.36 

43 

1.9 

41.1 

82 

0.71 

6.68 

6.38 

122 

22 

27 

0.78 

60 

8.9 

51.1 

125 

1.68 

13.03 

11.09 

123 

24 

24 

1.26 

97 

22.9 

74.1 

102 

2.60 

17.70 

13.52 

131 

22 

20 

0.34 

36 

1.7 

34.3 

93 

0.73 

7.82 

7.44 

132 

26 

20 

0.56 

54 

4.6 

49.4 

78 

1.11 

8.39 

7.67 

133 

24 

21 

1.00 

69 

14.5 

54.5 

89 

2.06 

12.59 

9.94 

201 

23 

40 

0.11 

29 

0.2 

28.8 

177 

0.23 

5.76 

5.72 

202 

25 

35 

0.80 

66 

9.4 

56.6 

143 

1.62 

11.10 

9.52 

203 

24 

40 

1.24 

92 

22.2 

69.8 

170 

2.56 

16.79 

12.74 

301 

36 

27 

0.12 

48 

0.2 

47.8 

76 

0.20 

3.89 

3.87 

302 

27 

29 

0.34 

47 

1.7 

45.3 

109 

0.66 

6.78 

6.53 

303 

32 

30 

0.57 

71 

4.8 

66.2 

96 

1.02 

7.29 

6.79 

311 

31 

25 

0.79 

76 

9.1 

66.9 

82 

1.43 

8.31 

7.31 

312 

34 

28 

1.05 

111 

16.0 

95.0 

84 

1.82 

10.09 

8.64 

313 

29 

27 

1.24 

109 

22.2 

86.8 

95 

2.32 

13.62 

10.85 

401 

40 

45 

0.12 

77 

0.2 

76.8 

115 

0.18 

5.06 

5.04 

402 

40 

45 

0.36 

96 

1.9 

94.1 

115 

0.57 

6.31 

6.18 

403 

40 

45 

0.52 

123 

4.0 

119.0 

115 

0.84 

8.08 

7.81 

411 

35 

60 

0.55 

90 

4.5 

85.5 

110 

1.31 

9.44 

8.74 

412 

37 

40 

0.79 

123 

9.1 

113.9 

142 

1.68 

10.87 

9.69 

413 

36 

50 

1.00 

134 

14.5 

119.5 

175 

0.94 

7.72 

7.34 

421 

42 

43 

0.11 

94 

0.2 

93.8 

104 

0.17 

5.60 

5.59 

422 

44 

56 

0.56 

144 

4.6 

139.4 

130 

0.85 

7.82 

7.57 

423 

42 

47 

0.99 

188 

14.2 

173.8 

114 

1.54 

11.20 

10.35 

analysis  techniques.  The  nondimensional  ice  resistance,  R,=  Rit/yBh*or  R./yBh*  is  consid¬ 
ered  to  be  primarily  a  function  of  Froude  number,  F  ,  based  on  ice  thickness  (Fn  =  V/'JgJi) 
and  of  the  nondimensional  ice  strength,  Cn  =  cr/y/r,  where  B  is  the  beam  at  the  waterline, 
y  is  the  specific  weight  of  water  and  h.  and  <x  are  the  thickness  and  flexural  strer  gth  of  ice, 
respectively.  Furthermore,  R,  is  assumed  to  be  a  linear  function  of  Cn  and  either  a  first-  or 
a  second-degree  polynomial  in  F  ,  that  is 

R,=  a  F2+  b  F  +  c  +  d  C  (4) 

n  n  n 

or 

R,  =  b  F  +c  +  dC  (5) 

Regression  analyses  in  accordance  with  eq  4  or  5  led  to  the  following  results.  For  the  total 
ice  resistance,  R  . 

R,=F2+  2.50  F  +2.68  +  0.0174  C  (r-r.9 61)  (4a) 

R,  =  5.21  F  +1.16  +  0.0195  C  (r  =  0.952)  (5a) 

*  n  n 

and  for  the  net  ice  resistance,  R 

'  1 
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R,(/yBh,  Measured  R,,/yBh, 


R.  =  0.22  F  2+  2.91  F  +  2.74  +  0.0157  C  (r  =  0.944) 

n  p  n 


(4b) 


R,  =  3.50  Fn  +  2.41  +  0.0162  Cn  (r  =  0.943).  (5b) 

where  r  is  the  correlation  coefficient.  The  correlation  coefficients  for  the  two  forms  of  the 
regression  equation  being  nearly  identical  indicates  that  these  two  forms  (eq  4  and  5)  are 
practically  indistinguishable.  The  data  are  shown  graphically  in  Figures  3a  and  3b  together 
with  the  regression  curves,  namely  eq  4a  for  the  total  resistance  and  eq  5b  for  the  net  re¬ 
sistance.  In  Figures  3c  and  3d  the  calculated  values  of  the  total  and  net  resistance,  re¬ 
spectively,  are  plotted  against  their  measured  values.  It  is  shown  that  all  measured  data  fall 
within  ±15%  of  the  calculated  values.  As  in  any  regression  analysis,  the  values  of  the 
coefficients  in  the  regression  equations  are  valid  only  for  the  test  range  of  the  variables,  that 
is,  for  the  present  resistance  tests 

0.17<  F  <2.6 


and 


60  <  C  <  180. 

n 


a.  Total  ice  resistance  vs  Froude  number.  b.  Net  ice  resistance  vs  Froude  number. 


c.  Meas •  r  .  -s  calculated  total  resistance.  d.  Measured  vs  calculated  net  resistance. 

Figure  3.  Results  of  resistance  tests  in  level  ice. 


Comparison  with  results  from  smooth  model 

The  results  of  the  resistance  tests  in  ice  with  the  smooth  R-class  model  (ice-hull  friction 
factor  of  0.04)  have  been  reported  earlier  (Tatinclaux  1 984).  These  data  were  reanalyzed  in 
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Rit/yBh,  Rlt/yBh 


terms  of  eq  4  and  5  with  the  following  results.  For  total  resistance  in  ice,  R 


R,  =  1.46  F2+  0.92  F  +  0.07  +  0.0376  C  (r  =  0.980)  (4c) 

w  n  n  it 

R,  =  4.75  F  - 1.48  +  0.0381  L  (r  =  0.964)  (5c) 

*  n  n 

and  for  net  ice  resistance,  R 

1 

R,  =  0.78  F2+  1.09  F  +  0.19  +  0.0364  C  (r  =  0.966)  (4d) 

R,  =  3.13  F  -  0.64  +  0.0367  C  (r  =  0.958).  (5d) 

*  n  n 


The  data  points  for  the  total  resistance  of  both  the  smooth  and  roughened  model  are 
plotted  in  Figure  4a,  while  those  for  the  net  resistance  are  plotted  in  Figure  4c.  Equations  4a 
and  4c  for  the  total  ice  resistance  of  both  models  are  compared  in  Figure  4b  for  the  two  values 
of  Cnof  60  and  180,  the  extremes  of  the  test  range  for  Cn-  Equations  4b  and  4d  for  the  net  ice 
resistance  are  compared  in  Figure  4d.  At  Cn  =  60,  the  roughened  model  exhibits  a  higher 
resistance  than  the  smooth  model,  as  expected,  but  at  C  =180  both  models  have  practically 
the  same  resistance.  This  finding  would  indicate  that  for  high  ice  strength,  or  thin  ice,  or 
both,  the  frictional  component  is  small  compared  to  the  breaking  component  of  the  total 
resistance.  Increasing  hull  roughness  has  then  relatively  little  effect  on  the  resistance. 
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o 
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1 
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a.  Total  ice  resistance  data. 
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Ik  Regression  function  for  total  resistance. 
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.  °.  — 

- 

o  .  cf. 

O  %  •  O  0  • 

•  _ 

o  °  O 

A 

• 

o 

-<b 

03 

J _ 1 _ 1 
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c.  Net  ice  resistance  data. 


if.  Regression  function  for  net  resistance. 


Figured.  Comparison  between  resistance  of  smooth  and  roughened  models  (CRREL  data). 
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PROPULSION  TESTS  WITH  ROUGHENED  MODEL 
Experimental  setup 

For  the  propulsion  tests,  the  R-class  model  was  equipped  with  two  thrust  and  torque 
dynamometers,  one  per  propeller  shatt,  rated  at  1100  N  in  thrust  and  11.3  Nm  in  torque. 
Both  propeller  shafts  were  driven  by  a  single  745-W  (1-hp)  variable-speed  electric  dc  motor 
that  was  equipped  with  a  tachometer  servo-mechanism  to  maintain  the  rotational  speed  at 
the  set  value.  A 1 : 1 .7  gear  reducer  between  the  output  shaft  of  the  motor  and  the  input  shafts 
of  the  dynamometers  limited  the  maximum  rotational  speed  of  the  propeller  shafts  to  about 
1000  rpm,  the  maximum  allowable  speed  for  the  dynamometers.  The  shaft  speed  was 
measured  by  a  magnetic  pickup  mounted  over  a  60-tooth  gear  fastened  to  the  portside 
propeller  shaft.  The  frequency  of  the  magnetic  pickup  was  converted  to  a  dc  voltage  bv  a 
frequency-to-voltage  converter  for  digital  sampling  by  the  data-acquisition  system. 

Test  procedure 

All  propulsion  tests  were  made  by  the  captive  model  technique,  whereby  the  ship  model 
was  connected  to  the  towing  carriage  by  the  same  tow-post  -force-block  arrangement  used 
in  resistance  tests.  Three  tests  were  made  in  each  ice  sheet  at  the  same  carriage  speed  but  for 
three  propeller  speeds.  In  each  test,  the  initial  values  of  the  thrust  and  torque  from  the 
dynamometers  and  the  initial  values  from  the  force  block  were  measured  with  the 
propellers  running  at  a  very  low  speed  (about  25  rpm)  and  with  the  ship  model  at  rest.  In 
this  way,  the  effect  of  friction  between  propeller  shafts  and  bushings  on  the  measurements 
of  thrust  and,  especially,  torque  wras  greatly  reduced  if  not  eliminated.  Also,  torque,  thrust 
and,  consequently,  pull  of  the  propellers  are  negligible  at  such  low  propeller  speed.  With 
the  ship  model  still  at  rest,  the  propeller  speed  was  increased  to  the  selected  value  for  the 
particular  test  and  the  data  acquisition  program  started. 

After  approximately  10  seconds  under  these  bollard  conditions,  the  carriage  was  set  in 
motion  at  the  selected  velocity  and  measurements  were  taken  over  about  one-third  of  the 
tank  length.  Once  the  carriage  had  stopped,  data  were  gathered  under  bollard  conditions 
for  an  additional  10  seconds,  the  propeller  speed  was  then  reduced  to  the  initial  low  value 
of  about  25  rpm  and  final  values  were  taken  for  comparison  with  the  initial  values, 
especially  values  for  thrust  and  torque.  The  ship  model  was  backed  into  the  broken  channel 
and  the  procedure  repeated  for  the  next  test  at  the  same  carriage  speed  and  a  different 
propeller  speed.  The  backing  distance  was  such  that  in  the  next  test  the  ship  model  reached 
steady  speed  before  entering  the  unbroken  ice  sheet. 

The  purpose  of  the  bollard  tests  at  the  beginning  and  the  end  of  each  propulsion  test  was 
to  check  the  proper  functioning  of  the  dynamometers.  Throughout  the  tests,  the  analog 
signals  from  the  force  block,  the  thrust  and  torque  outputs  of  one  dynamometer,  and  the 
propeller  speed  counter  were  monitored  on  a  four-channel  chart  recorder  to  ensure  visually 
that  all  systems  were  functioning  properly. 

In  each  series  of  three  tests  per  ice  sheet,  the  propeller  speed  for  the  first  test  was  selected 
so  that  the  ship  model  would  be  underpropelled,  that  is  the  pull  would  be  negative;  in  the 
second  test,  the  propeller  speed  was  increased  sufficiently  to  ensure  that  the  model  would 
be  overpropelled  (positive  pull);  and  in  the  third  test  the  propeller  speed  was  set  at  an 
intermediate  value  in  an  attempt  to  have  a  pull  nearly  equal  to  zero. 

In  addition  to  the  tests  in  ice,  propulsion  tests  in  clear  water  were  also  made  for 
comparison  with  the  test  results  reported  by  Murdey  (1980).  In  these  tests,  the  model  was 
towed  at  constant  velocity  and  the  propeller  speed  was  increased  at  regular  intervals  along 
the  length  of  the  tank. 

The  dynamometers  and  the  force  block  were  calibrated  by  conducting  bollard  tests  in 
clear  water  over  a  wide  range  of  propeller  speeds  and  adjusting  the  respective  calibration 
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<N)  RPM  Torque  (Nm)  Ve  I  Com/*) 


coefficients  to  match  the  propeller  performance  reported  by  Murdev  (1980).  The  oollnrd 
thrust  and  torque  for  each  propeller  and  the  total  pull  were  found  to  be  proportional  to  the 
square  of  the  propeller  speed,  namely,  when  n  is  expressed  in  revolutions  per  second 


Pull  (N)  =  1.175  id 

(6) 

T  (N)  =  0.61 

a  a 

(7) 

Q^Nm)  =0.0199  iO 

(8) 

which  correspond  to  thrust  and  torque  coefficients  of  K  =0.339  and  K  =0.0402,  and  a  thrust 
deduction  factor  of  t  =  0.04. 

Data  acquisition  system 

Data  were  acquired  on  a  NEFF  620  system,  consisting  of  a  Model  100  signal  conditioner 
and  a  Model  300  signal  processor  (analog-to-digital  converter),  controlled  bv  a  HP-9845B 
desktop  computer.  Seven  data  channels  were  scanned  and  sampled,  namely  carriage  speed, 
force  block  (pull),  propeller  speed,  and  thrust  and  torque  of  both  propellers.  The  sampling 
rate,  S  (in  milliseconds),  of  theanalog  signals  was  selected  so  that  a  minimum  of  32  samples 
were  taken  per  propeller  revolution  (8  samples  per  propeller  blade),  that  is 

2<S(ms)<  1000/(32;/,). 

The  data  in  digitized  form  were  stored  on  floppy  disks  for  subsequent  analysis.  An 
example  of  the  data  signals  is  shown  in  Figure  5. 


a.  Carriage  velocity,  starboard  torque,  propeller  speed  mid  pull. 


Figure  5.  Example  of  data  signals  for  test  no.  $10  (h  =  3.5  nun;  c  =  20  kPa,  V  =  0.35  m/s;  n  t  =  377  ;/>m). 
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b.  Thrust  and  torque  on  both  propellers. 

Figure  5  (cant'd). 

Data  presentation  and  analysis 

Data  presentation 

For  each  test,  the  data  were  averaged  over  the  period  of  uniform  carriage  velocity.  The 
test  data  in  dimensional  form  are  listed  in  Table  5  and  in  nondimensional  form  in  Table  6, 
where  /g  is  the  appparent  advance  coefficient,  K(  the  thrust  coefficient  and  the  torque 
coefficient,  defined  as 

l=V/nD  (9) 

K  =  T  /(2p  »2  D4)  (10) 

t  a  1  3 

=  Qa/ (2  pnj  D5)  (11) 

where  p  is  water  density  and  D  propeller  diameter. 

Self-propulsion  points 

The  results  of  the  three  propulsion  tests  carried  out  at  a  given  carriage  velocity  in  one  ice 
sheet  were  interpolated  to  the  point  of  zero  pull  to  obtain  the  self-propulsion  point  for  the 
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Table  5.  Results  of  propulsion  tests  with  roughened  R-class 
model. 


Pari  Starboard 


Test 

tto. 

h 

t 

(mm) 

a 

(kh) 

V 

(ml s) 

n 

a 

(rpnt) 

Pull 

(N) 

T 

n 

(N) 

<2* 

(Km) 

T 

a 

(N) 

Q„ 

(Nm) 

a.  Tests  iu  level  ice 

511 

38 

45 

0.34 

398 

-72 

21.7 

0.78 

21.5 

0.87 

512 

600 

-11 

55.2 

1.42 

53.6 

1.47 

513 

653 

21 

64.9 

1.64 

64.6 

1.64 

552 

31 

45 

0.52 

400 

-47 

17.9 

0.75 

17.9 

0.61 

551 

500 

NA 

32.9 

1.05 

27.8 

1.51 

553 

702 

41 

67.2 

1.93 

68.3 

1.98 

521 

37 

32 

0.74 

599 

-47 

38.1 

1.97 

42.9 

1.19 

522 

708 

-23 

55.4 

2.34 

59.1 

2.16 

523 

757 

10 

69.9 

2.34 

69.9 

2.37 

543 

25 

25 

0.52 

393 

-15 

18.0 

0.62 

17.9 

0.62 

541 

497 

0 

32.3 

0.98 

30.7 

1.04 

542 

602 

37 

49.8 

1.37 

52.1 

1.34 

531 

26 

26 

0.76 

498 

-35 

22.9 

1.00 

21.8 

1.21 

532 

0.75 

600 

-5 

35.6 

1.53 

38.9 

1.56 

533 

0.74 

625 

8 

45.3 

1.52 

45.6 

1.56 

561 

22 

33 

1.24 

707 

-27 

39.5 

1.74 

41.5 

1.68 

562 

1.20 

760 

-10 

54.1 

2.32 

54.0 

2.39 

653 

1.20 

812 

6 

64.1 

2.47 

66.5 

2.45 

b.  Tests  in  clear  water 

911 

0.34 

541 

78 

43.0 

1.01 

42.6 

1.02 

912 

482 

60 

33.7 

0.81 

34.1 

0.82 

913 

410 

42 

25.4 

0.56 

23.7 

0.55 

914 

343 

25 

14.9 

0.36 

15.8 

0.37 

915 

275 

16 

9.1 

0.21 

9.5 

0.23 

916 

211 

5 

4.7 

0.12 

4.8 

0.09 

921 

0.52 

541 

68 

40.8 

1.00 

37.1 

0.93 

922 

483 

50 

30.1 

0.77 

28.9 

0.71 

923 

418 

39 

21.3 

0.55 

20.9 

0.51 

924 

338 

21 

12.6 

0.31 

12.6 

0.30 

925 

262 

8 

6.1 

0.16 

6.6 

0.17 

931 

0.74 

707 

110 

64.5 

1.67 

61.2 

1.57 

932 

604 

69 

45.4 

1.19 

43.9 

1.11 

933 

515 

43 

29.5 

0.81 

29.4 

0.75 

934 

410 

20 

15.0 

0.43 

16.1 

0.39 

935 

325 

6 

6.5 

0.22 

7.6 

0.21 

941 

1.20 

749 

59 

52.2 

1.52 

55.0 

1.55 

942 

658 

39 

39.4 

1.13 

40.4 

1.14 

943 

555 

13 

24.0 

0.63 

25.0 

0.68 

12 


Table  6.  Nondimensional  form  of 
propulsion  test  results. 


Test  no. 

F 

n 

C 

It 

K, 

a.  TVs 

Is  in  te 

I'd  ice 

511 

0.56 

121 

0.25 

0.274 

0.0505 

512 

0.17 

0.302 

0.0390 

513 

0.15 

0.304 

0.0373 

552 

0.94 

148 

0.3S 

0.224 

0.0412 

553 

0.30 

0.243 

0.0497 

551 

0.22 

0.275 

0.0385 

521 

1.23 

88 

0.36 

0.226 

0.0427 

522 

0.30 

0.228 

0.0436 

523 

0.28 

0.244 

0.0399 

543 

1.05 

102 

0.39 

0.232 

0.0390 

541 

0.30 

0.255 

0.0397 

542 

0.25 

0.281 

0.0363 

53! 

1.50 

102 

0.44 

0.180 

0.0432 

532 

1.49 

0.36 

0.207 

0.0416 

533 

1.47 

0.34 

0.233 

0.0383 

561 

2.67 

153 

0.51 

0.162 

0.0332 

562 

2.58 

0.46 

0.187 

0.0396 

563 

2.58 

0.43 

0.198 

0.0362 

h  Tt 

•>/s  in , 

clear  ?  enter 

911 

0.18 

0.292 

0.0337 

912 

0.21 

0.292 

0.0340 

913 

0.24 

0.292 

0.0320 

914 

0.29 

0.261 

0.0301 

915 

0.36 

0.246 

0.0282 

916 

0.47 

0.213 

0.0229 

921 

0.28 

0.266 

0.0320 

922 

0.32 

0.253 

0.0308 

923 

0.36 

0.241 

0.0294 

924 

0.45 

0.221 

0.0259 

925 

0.58 

0.185 

0.0233 

931 

0.31 

0.251 

0.0315 

932 

0.36 

0.245 

0.0306 

933 

0.42 

0.222 

0.0285 

934 

0.53 

0.185 

0.0237 

935 

0.66 

0.133 

0.0198 

941 

0.47 

0.191 

0.0266 

942 

0.53 

0.184 

0.0254 

943 

0.63 

0.159 

0.0206 

O  50  100  150 


r  p  m 

Figure  6.  Example  of  determination  of  self-propulsion  point 
by  interpolation  of  propulsion  test  results  (test  series  510). 


n  4 , 


Study 

Overload  Tests  (Murde<  '9Q0  ’ 
Sei’-propulS'Of  •  C>eor  /.3*e’  iV/ie. 


Figure  7.  Thrust  deduction  t actors  from  several 
studies. 
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Table  7,  Self-propulsion  points  in  ice  of  roughened  R-class  model. 


Test  h.  0i  V  R(.(  na  Total  Total 

series  (mm)  (kPa)  (ml s)  (N)  (rpm)  T/N)  QJNm)  Ja  K(  10  t 


510 

38 

45 

0.34 

93 

620 

116 

3.01 

0.16 

0.304 

0382 

020 

550 

31 

45 

052 

84 

643 

111 

3.43 

0.24 

0269 

0.404 

0.24 

520 

37 

32 

0.74 

117 

744 

133 

4.71 

0.29 

0240 

0.413 

0.12 

540 

25 

25 

0.52 

50 

497 

65 

2.02 

0.30 

0.262 

0.397 

0.23 

530 

26 

26 

0.75 

70 

610 

82 

3.05 

0.36 

0.210 

0.395 

0.15 

560 

22 

33 

1.20 

92 

786 

122 

4.78 

0.45 

0.197 

0.375 

0.25 

Note:  Total  resistance  Rj(=  Ri  +  Rqw,  with  R.  given  by  eq  4b  and  RQwby  eq  3. 


corresponding  ice  thickness,  ice  strength  and  ship  speed.  An  example  of  this  interpolation 
procedure  is  shown  in  Figure  6.  The  self-propulsion  points  thus  obtained  are  given  in  Table 
7,  where  R.  =  R  +  Rowl  Rt was  calculated  according  to  eq4b  and  the  thrust  deduction  factor 
was  defined  by 

f  =  l-J?it/T.  (12) 

The  accuracy  of  the  thrust  deduction  factor,  t,  depends  on  the  accuracy  of  the  resistance 
estimate  and  of  the  thrust  measurements.  For  given  ice  conditions  of  thickness  and  flexural 
strength,  ship  resistance  in  level  ice,  in  contrast  with  that  in  ice-free  water,  can  usually  be 
estimated  at  best  within  ±1 0%  of  its  actual  value.  Therefore,  even  if  the  propeller  thrusts  are 
accurately  measured,  significant  uncertainties  remain  in  the  values  of  t  calculated  from  the 
results  of  independent  resistance  and  propulsion  tests.  Figure  7  presents  the  values  of  t  at 
the  self-propulsion  points  as  determined  in  the  present  study.  Also  shown  in  Figure  7  are 
the  values  of  t  evaluated  from  the  results  of  overload  tests  in  ice-free  water  (Murdey  1980) 
for  the  conditions  of  model  speed  and  propeller  speed  corresponding  to  the  self-propulsion 
points  in  ice,  and  the  values  of  t  obtained  by  Murdey  (1980)  at  the  self-propulsion  points  in 
clear  water.  As  can  be  seen  from  Figure  7,  the  values  of  t  in  level  ice  are  scattered  about  those 
in  ice-free  water. 

Until  a  better  method  to  evaluate  the  thrust  deduction  factor  in  ice  is  available,  it  is 
recommended  that  the  values  obtained  under  clear  water  conditions  be  used.  In  the  present 
case,  they  can  be  expressed  by 

t  =  0.04  +  0.12  V  (13) 

which  is  valid  for  V  <  1 .20  m/ s,  the  maximum  speed  in  the  present  test  series.  If  it  is  assumed 
that  the  thrust  deduction  is  the  same  at  full  scale  and  at  model  scale  for  corresponding  ship 
and  model  speeds,  then  the  expression  for  t  at  full  scale  is 

t  =  0.04  +  0.027  V.  (14) 

Equation  14  is  valid  for  V  <  5.37  m/s  at  full  scale. 

The  results  of  propulsion  tests  with  the  smooth  model  had  led  to  the  conclusion  that  the 
thrust  deduction  factor,  f,  in  level  ice  was  constant,  independent  of  velocity  and  equal  to 
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about  0.20.  Because  of  doubts  about  the  proper  functioning  of  the  thrust  dynamometers 
during  this  previous  study,  as  explained  below,  eq  13  and  14  are  assumed  to  apply  also  to 
the  results  with  the  smooth  model  in  the  following  analyses  and  extrapolation  to  full-scale 
conditions. 

Propeller  coefficients 

The  thrust  and  torque  coefficients,  K  and  K  ,  for  all  propulsion  tests  both  in  level  ice  and 
in  ice-free  water  are  plotted  versus  /a  on  Figure  8,  where  the  values  of  Kf  and  K  obtained 
with  the  smooth  model  are  also  shown.  It  can  be  seen  that  the  thrust  coefficient  for  the 
rougher  model  in  level  ice  is  only  slightly  smaller  than  that  in  ice-free  water.  K(  for  the 
smooth  model  from  the  previous  study  was  generally  higher  than  for  both  the  rougher 
model  and  the  ice-free  conditions.  This  latter  observation  would  indicate  inaccuracies  in  or 
malfunction  of  the  thrust  dynamometer  during  the  previous  study.  Therefore,  in  the 
following  analyses  and  full-scale  extrapolations,  the  thrust  coefficients  obtained  from  the 
present  propulsion  tests  were  assumed  valid  for  both  the  rough  and  smooth  model.  On  the 
other  hand,  the  values  of  ^obtained  in  both  studies  from  propulsion  tests  in  level  ice  are 
comparable  and  consistently  greater  than  those  for  the  ice-free  conditions.  The  changes  in 
K(and  ^between  level  ice  and  ice-free  conditions  are  attributed  to  ingestion  of  ice  floes  by 
the  propellers.  The  resulting  impacts  of  the  propeller  blades  on  ice  floes  increase  the  torque 
on  the  propellers  and  decrease  the  thrust  delivered  by  the  propellers.  The  experimental  data 
were  analyzed  using  linear  regression  and  forcing  the  regression  curves  to  pass  through  the 
respective  values  of  K  and  K  at  bollard.  For  ice-free  conditions,  the  resulting  equations 
were  q 


K=  0.339  -  0.303  /  (r  =  0.986) 

I  3 


(15) 


Kq  =  0.0402  -  0.0308  /a  (r  =  0.969).  (16) 

For  the  case  of  level  ice  conditions,  the  equations  were 

Kt  -  0.339  -  0.318  /a  (r  =  0.972)  (17) 

and 


K  =0.0402-0.0066/  (r  =  0.38). 

q  a 


These  equations  are  shown  in  Figure  8. 


a.  Thrust  coefficient. 


(18) 


b.  Torque  coefficient. 


Figure  8.  Propulsion  coefficients  vs  apparent  advance  coefficient. 
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FULL-SCALE  PERFORMANCE  PREDICTIONS 
AND  COMPARISON  WITH  EXISTING  DATA 


Performance  predictions 

Equations  5b,  4d  and  14  to  1 8  are  the  basis  for  predicting  full-scale  performance  from  the 
model  test  results.  For  the  hull  and  propeller  characteristics  of  the  R-class  icebreaker,  these 
equations  become 


Rj  (kN)  =  212.3  Vh\5  +458.0 It] 

+  0.314  o,  h  j  (rough) 

(19a) 

R;  (kN)  =  15.1  V2/i;  +66.1  VU\ 

3  +  36.1  l;2  +  0.705  Oj  h  j  (smooth) 

(19b) 

R.  (kN)  =  R.  +  R 

It  1  OW 

(20) 

T  (kN)  =  Rit/(0.96  -  0.027  V) 

(21) 

T  (kN)  =  195  it2- 44.5  V  n 

a  '  a  a 

(level  ice) 

(22a) 

T  (kN)  =  195  n2- 42.4  V  n 

a  a  a 

(clear  water) 

(22b) 

Qa(kNm)  =  95.4  n\-  3.80  1/  »a 

(level  ice) 

(23a) 

Qa(kNm)  =  95.4  n\-  17.8  V  »a 

(clear  water) 

(23b) 

PD  (kW)  -2k  n aQa 

(24) 

In  eq20,  R  is  the  full-scale  resistance  in  ice-free  water.  This  resistance  was  predicted  from 
model  test  experiments  by  Murdey  (1980).  For  the  purpose  of  comparing  the  predicted  ship 
performance  in  level  ice  with  existing  full-scale  trial  data,  Murdey's  results  were  expressed 
in  the  following  form 


(kN)  =  4.41  V  +  2.56  V 2  +  0.155  V  3 

(V  <6  m/s) 

(25a) 

(kN)  =  152  +  54  (V-6)  +  14.2  (V-6)3 

(V  >  6  m/s) 

(25b) 

In  eq  19-25,  h.  is  in  meters,  V  in  meters  per  second,  o.in  kilopascals,  and  n  in  s"1  (revolutions 
per  second).  It  should  be  reiterated  that  these  equations  are  strictly  valid  within  the 
experimental  range 0.17-2.6  for  the  Froude  number,  F  ,  and  within  the  range  60-180  for  the 
dimensionless  ice  strength,  C  . 

Predicted  total  resistance,  total  thrust  and  shaft  power  for  the  conditions  It.  -  0.45  and 
0.70  m  and  a  =  450  and  600  kPa  are  presented  in  Table  8  for  both  models  and  plotted  in 
Figure  9  for  the  smooth  model  and  in  Figure  1 0  for  the  rough  model.  Conversely,  given  shaft 
power  and  ice  strength,  eq  19-25  can  be  used  to  predict  the  speed  at  which  the  vessel  will 
break  level  ice  of  given  thickness.  Such  performance  predictions  of  V  vs  It.  at  full  power  of 
11,000  kW  for  ice  strength  of  450  kPa  are  listed,  with  corresponding  ice  resistance  and  total 
thrust,  in  Table  9a  where  the  values  of  K(and  K  from  propulsion  tests  in  level  ice  were  used, 
and  in  Table  9b  where  those  from  tests  in  clear  water  were  used.  The  tables  contain 
predictions  based  on  the  resistance  equations  derived  from  the  resistance  test  results  with 
the  smooth  model  and  the  roughened  model.  These  predicted  performances  are  shown  as 
maximum  velocity  at  full  power  versus  ice  thickness  in  Figure  1 1  a  for  the  smooth  model  and 
in  Figure  lib  for  the  roughened  model. 


( MN)  ,  P/ 10  (MW)  R.  T  (MN|  P/IOIMW] 


Table  8.  Predicted  performance  of  both  models  for  given  ice  condi¬ 
tions  (K,  and  K  from  model  tests  in  ice). 

'  q 


Ice  Smooth  model  _ Rough  model 


conditions 

\(kn) 

RJMN) 

T  (MN) 

a 

PDf  MW) 

R  ,(MN) 

it 

T  (MN) 

a 

PD(MW) 

3 

0.21 

0.23 

1.18 

0.27 

0.29 

1.61 

4 

0.24 

0.27 

1.64 

0.31 

0.34 

2.23 

h  =0.45  m 

5 

0.28 

0.31 

2.25 

0.35 

0.40 

3.00 

a  =450  kPa 

6 

0.32 

0.36 

3.04 

0.40 

0.45 

3.93 

7 

0.37 

0.42 

4.05 

0.44 

0.51 

5.04 

8 

0.42 

0.49 

5.31 

0.49 

0.58 

6.36 

3 

0.34 

0.37 

2.18 

0.53 

0.58 

3.99 

4 

0.39 

0.43 

2.95 

0.60 

0.66 

5.27 

h  =070  m 

5 

0.44 

0.49 

3.95 

0.67 

0.76 

6.81 

a  =450  kPa 

6 

0.50 

0.57 

5.22 

0.75 

0.86 

8.61 

7 

0.57 

0.66 

6.81 

0.83 

0.96 

10.70 

8 

0.65 

0.77 

8.78 

0.91 

1.07 

13.12 

3 

0.26 

0.28 

1.53 

0.29 

0.32 

1.78 

4 

0.29 

0.32 

2.04 

0.33 

0.37 

2.43 

h  =0.45  m 

5 

0.32 

0.36 

2.71 

0.37 

0.42 

3.22 

o  =600  kPa 

6 

0.37 

0.42 

3.57 

0.42 

0.48 

4.18 

7 

0.41 

0.48 

4.64 

0.46 

0.54 

5.32 

8 

0.47 

0.55 

5.99 

0.51 

0.61 

6.67 

3 

0.41 

0.45 

2.85 

0.56 

0.61 

4.33 

4 

0.46 

0.51 

3.71 

0.63 

0.70 

5.66 

h  =0.70  m 

5 

0.51 

0.58 

4.81 

0.71 

0.79 

7.24 

a  =600  kPa 

6 

0.58 

0.66 

6.19 

0.78 

0.89 

9.09 

7 

0.65 

0.75 

7.89 

0.86 

1.00 

11.24 

8 

0.72 

0.85 

9.98 

0.94 

1.11 

13.71 

V  (m  /$  )  V  ( m  /s  ) 

c.  h  =  0.45  in;  o  =  600  kPa.  d.  h  =  0.70  m;  a  =  600  kPa. 

I  I  II 

Figure  9.  Predicted  performance  in  lei'el  ice  from  test  results  with  smooth  model  (Kpind  K  from  tests  in  level  ice). 
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Table  9.  Predicted  performance  in  level  ice  at  full  power. 


Smooth  model 


Roiofh  model 


O. 

h 

1 

R„ 

T 

it 

R, 

T( 

(kPa) 

(in) 

(kn) 

(MN) 

(MN) 

( km 

(MN) 

(MN) 

a.  K fiimi  K^  from  model  te* 

■Is  in  ice. 

450 

0.20 

13.7 

0.48 

0.60 

14.4 

0.45 

0.56 

0.30 

12.6 

0.54 

0.67 

13.0 

0.51 

0.64 

0.40 

11.4 

0.59 

0.73 

11.4 

0.59 

0.74 

0.50 

10.5 

0.64 

0.78 

9.9 

0.68 

0.82 

0.60 

9.7 

0.69 

0.83 

8.4 

0.76 

0.90 

0.70 

8.9 

0.73 

0.87 

7.1 

0.84 

0.97 

0.80 

8.2 

0.77 

0.91 

5.9 

0.91 

1.04 

0.90 

7.6 

0.81 

0.95 

4.9 

0.98 

1.10 

1.00 

7.0 

0.84 

0.98 

3.9 

1.04 

1.15 

600 

0.20 

13.5 

0.49 

0.61 

14.3 

0.46 

0.57 

0.30 

12.3 

0.55 

0.69 

12.9 

0.52 

0.65 

0.40 

11.1 

0.61 

0.75 

11.2 

0.60 

0.74 

0.50 

10.1 

0.66 

0.81 

9.7 

0.69 

0.83 

0.60 

9.2 

0.71 

0.86 

8.2 

0.77 

0.9 1 

0.70 

8.4 

0.76 

0.90 

6.9 

0.85 

0.99 

0.80 

7.7 

0.80 

0.94 

5.7 

0.93 

1.05 

0.90 

7.0 

0.84 

0.98 

4.6 

1 .00 

111 

1.00 

6.4 

0.88 

1.01 

3.6 

1 .06 

1.17 

/>. 

K t,nht  K 

from  moiU'l  tests  in 

clear  woler. 

450 

0.20 

16.0 

0.75 

0.94 

16.6 

0.74 

0.92 

0.30 

15.0 

0.77 

0.97 

15.7 

0.76 

0.94 

0.40 

13.9 

0.80 

0.99 

14.3 

0.79 

0.98 

0.50 

12.7 

0.82 

1.03 

12.3 

0.83 

1.04 

0.60 

11.5 

0.85 

1.06 

10.3 

0.89 

1 .09 

0.70 

10.6 

0.88 

1 .09 

8.6 

0.96 

1.14 

0.80 

9.7 

0.92 

1.11 

7.1 

1.02 

1.18 

0.90 

8.9 

0.95 

1.13 

5.7 

1.07 

1.22 

1.00 

8.2 

0.97 

1.15 

4.5 

1.12 

1.25 

6(H) 

0.20 

15.9 

0.75 

0.94 

1 6.6 

0.74 

0.92 

0.30 

14.8 

0.77 

0.97 

15.6 

0.76 

0.95 

0.40 

13.6 

0.80 

1.00 

14.2 

0.79 

0.99 

0.50 

12.2 

(V83 

1.04 

12.1 

0.84 

1.05 

0.60 

11.0 

0.87 

1 .08 

10.1 

0.90 

1.10 

0.70 

10.0 

0.9 1 

1.10 

8.4 

0.97 

1.15 

0.80 

9.1 

0.94 

1.13 

6.8 

1.03 

1.19 

0.90 

S3 

0.97 

1.15 

5.4 

1.08 

1.22 

1 .00 

73 

1.00 

1.17 

4.2 

1.13 

1.26 

6 


a.  h  =  0.45  in ;  a  =  450  kPn.  b.  h  =  0.71)  in;  a  =  450  kPn. 

Figure  10.  Predicted  performance  in  level  ieefwiu  test  results  with  roughened  model  (K  pmd  K  from  tests  in  level 
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{  M  N  )  ,  P/IOIMW) 


c.  h  =  0.45  m;  a  =  600  kPa. 

I  I 

Figure  10  (coat'd). 


d.  h  =  0.70  m;  a  =  600  kPa. 

1  l 


a.  Smooth  model.  b.  Roughened  model. 


Figure  1 1 .  Predicted  performance:  V  i's  h  at  constant  power  (with  a  =  450  kPa,  t  from  clear  water 
tests  and  K(  and  K  from  tests  in  level  ice  and  in  clear  -water). 

Comparison  between  predictions  and  measurements 

In  its  report  to  the  18th  International  Towing  Tank  Conference  (ITTC  1987),  the  Ice 
Committee  presented  a  reanalysis  of  available  full-scale  data  from  trials  with  the  CCGS 
Radisson  ( Edwards  et  al.  1 981 )  and  the  CCGS  Franklin  (Michailidis  and  Murdev  1 981 ).  These 
adjusted  full-scale  data,  corresponding  to  ice  strength  of  either  400  or  480  kPa,  are  given  in 
Table  10.  In  addition,  the  particular  set  of  full-scale  data  by  Michailidis  and  Murdey  (1981) 
are  listed  in  Table  1 1  as  they  appear  to  the  writers  as  being  the  best  available  set  of  field  data 
with  the  R-class  icebreaker.  Both  sets  of  data  are  plotted  in  Figures  11a  and  lib  for 
comparison  with  the  above  ship  performance  predictions. 

From  these  figures  it  can  be  seen  that  the  performance  predictions  based  on  the 
roughened  model  resistance  and  the  values  of  and  in  clear  water  are  in  very  good 
agreement  w.  h  the  full-scale  data,  especially  for  ice  thicknesses  less  than  0.8  m.  This 
observation  must  be  interpreted  as  a  strong  indication  that  effect  on  propeller  performance 
of  ice  ingestion  into  the  propeller  disks  is  far  greater  at  model  scale  than  at  full  scale.  This 
exaggerated  ice  effect  at  model  scale  can  be  caused  by  more,  and  relatively  larger,  ice  floes 
being  ingested  by  the  propellers  than  it  full  scale  because  of  incorrect  modeling  of  ice 
properties,  such  as  density,  bulk  modulus  and  fracture  toughness. 
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Table  10.  Full-scale  data  for  R- 
class  performance  in  level  ice 
(ITTC  1987). (Full-scale  ice  flex¬ 
ural  strength  was  407  or  480  kPa.) 

a.  Speed  in  level  ice  at  full 
power  <11, 000  kW). 


h  (iii)  V(kn) 

I 


0.2 

15.9 

0.3 

14.7 

0.4 

13.3 

0.5 

11.7 

0.6 

10.0 

0.7 

8.2 

0.8 

6.6 

0.9 

4.7 

1.0 

2.9 

b.  Total 

resifihmct 

■,  thrust  and  power  in  0.7- 

in  level 

lev. 

R. 

T  PD 

(km 

(MN) 

(MN)  (MWi 

3 

0.50 

0.54  3.43 

4 

0.54 

0.58  4.00 

5 

0.59 

0.65  4.93 

6 

0.65 

0.74  6.21 

7 

0.74 

0.85  7.75 

8 

0.85 

0.97  9.36 

Table  11,  Full-scale  trial  data  (after  Michailidis 
and  Murdey  1981). 


V(ku) 

h  (m) 

l 

h(M> 

a(kPa) 

rpm 

P(AIV1 

T  IkN) 

a 

Port  propeller 

2.18 

0.704 

0.08 

510 

115.7 

1840 

289.0 

4.64 

0.620 

0.13 

480 

126.3 

2320 

318.8 

5.38 

0.674 

0.08 

460 

132.6 

2560 

348.7 

4.05 

0.676 

0.15 

NA 

142.6 

3700 

413.5 

7.06 

0.650 

0.11 

NA 

155.7 

4180 

443.4 

8.23 

0.666 

0.12 

NA 

170.3 

5420 

528.1 

7.64 

0.613 

0.09 

NA 

162.4 

4610 

478.3 

Starboard  propeller 

2.18 

0.704 

0.08 

510 

122.0 

2020 

289.0 

4.64 

0.620 

0.13 

480 

127.0 

2190 

328.8 

5.38 

0.674 

0.08 

460 

133.9 

2950 

358.7 

4.05 

0.676 

0.15 

NA 

146.4 

3750 

413.5 

7.06 

0.650 

0.11 

NA 

154.7 

3840 

473.3 

8.23 

0.666 

0.12 

NA 

173.8 

5480 

523.1 

7.64 

0.613 

0.09 

NA 

164.6 

4580 

478.3 

COMPARISON  WITH  RESULTS  FROM  OTHER  FACILITIES 

As  mentioned  in  the  Introduction,  this  series  of  ship  model  tests  was  part  of  an  interna¬ 
tional  cooperative  effort  under  the  aegis  of  the  ITTC.  The  results  of  the  resistance  and 
propulsion  tests  made  with  the  smooth  model  by  the  various  participating  organizations 
were  presented  at  the  1 984  ITTC  in  Gdteborg,  Sweden  (ITTC  84) .  The  results  of  the  resistance 
tests  with  the  roughened  model  were  reported  by  the  Ice  Committee  to  the  18th  ITTC  in 
Kobe,  Japan  (ITTC  87).  The  results  of  the  propulsion  tests  with  the  latter  model,  which  have 
not  yet  been  completed  by  all  facilities,  are  to  be  presented  at  the  19th  ITTC  in  1990  in 
Madrid,  Spain.  The  results  of  the  resistance  and  propulsion  tests  on  both  the  smooth  and 
roughened  1:20  scale  R-class  model  by  the  various  facilities  are  listed  in  Appendix  A. 

Data  on  resistance  and  propulsion  performance  from  ship  model  tests  in  ice  are  analyzed 
by  using  standard  linear,  multilinear  or  nonlinear  regression  analysis  methods.  Since  there 
exists  no  satisfactory  analytical  formulation  of  the  interaction  between  ship  and  ice  to 
suggest  the  form  of  the  regression  equation,  the  choice  of  this  equation  is  arbitrary  and  left 
to  the  individual  researcher.  Based  on  their  experience,  the  writers  have  assumed  that  the 
selected  dependent  variable  (e  g.,  nondimensional  resistance,  thrust  or  torque  coefficient) 
could  be  expressed  as  a  first  or  second  degree  polynomial  of  the  independent  variables  (e.g., 
Froude  number,  nondimensional  ice  strength  or  advance  coefficient).  For  any  regression 
analysis  of  experimental  data  to  yield  reasonably  reliable  results,  the  data  set  should  be 
large  enough  and  more  than  one  test  should  have  been  made  for  the  same  experimental 
conditions  to  check  the  repeatablity  of  the  data.  Finally,  when  comparing  several  data  sets, 
the  range  of  the  independent  variables  should  be  nearly  the  same  from  one  data  set  to  an¬ 
other. 


20 


Table  12.  Range  of  resistance  test  conditions  with  smooth 


model 

Facility 

No.  of 
Tests 

Range  of  test  conditions 

h  (mm) 

t 

a(kPa) 

1 

Vtin/s) 

F 

n 

C 

II 

CRREL 

20 

20— tO 

25-55 

0.11-1.31 

0.19-2.85 

64-261 

HSVA 

25 

22-51 

18-43 

0.10-0.71 

0.16-1.45 

50-182 

JSR1 

33 

10-37 

27-65 

0.11-1.23 

0.19-3.95 

100-570 

NKK 

6 

35-36 

39-44 

0.11-1.24 

0.19-2.12 

114-125 

WARC 

7 

19-36 

19-32 

0.58-0.59 

0.98-1.34 

74-107 

Note:  Ice  thicknesses  were  rounded  to  the  nearest  millimeter,  and  ice 
strengths  to  the  nearest  kilopascal. 


Tests  with  smooth  model 

Resistance  tests 

Five  facilities,  CRREL,  HSV  A,  JSR1,  NKK  and  W  ARC,  conducted  resistance  tests  in  level 
ice  with  the  smooth  1:20  scale  R-class  model.  The  range  of  test  conditions  at  these  various 
facilities  is  listed  in  Table  12,  and  the  test  data  are  plotted  in  Figure  12. 

Note  that  both  NKK  and  WARC  conducted  only  a  very  limited  number  of  tests. 
Furthermore,  all  NKK  tests  were  made  at  practically  constant  ice  thickness  and  ice  strength 
(ice  strength  can  usually  be  measured  with  no  better  than  a  1 0%  accuracy),  and  velocity  was 
the  only  independent  variable.  On  the  contrary,  all  WARC  tests  were  made  at  a  single 
velocity,  and  at  only  two  ice  thicknesses  and  two  ice  strengths;  therefore,  the  corresponding 
data  are  insufficient  for  a  reliable  regress,  on  analysis  involving  more  than  one  independent 
variable. 

Statistical  analysis  of  all  test  data  indicated  that  the  model  total  ice  resistance,  Rjt,  was 
mainly  influenced  by  model  speed,  V,  and  ice  thickness,  It,  and  only  marginally  affected  by 
ice  strength,  a..  Consequently,  the  nondimensional  ice  resistance,  R,=  R  /yBIr,  is  primarily 
a  function  of  the  Froude  number,  Fn  =  V  / 'Igh  and  varies  only  weakly  with  C  =  ojyh.,  as 
is  evident  from  Figure  13a.  Table  13  shows  the  distribution  of  the  tests  with  the  smooth  R- 
class  model  among  various  Froude  number  ranges. 

Only  four  tests  were  conducted  at  Fn greater  than  2.75,  one  by  CRREL  and  three  by  JSR1. 
As  can  be  seen  from  Table  A1  and  Figure  13a,  each  of  these  four  tests  was  made  at  a  distinct 
value  of  Fn;  that  is,  the  repeatability  of  the  test  results  at  very  high  Fncannot  be  assured.  Since 
such  isolated  points  at  the  extreme  end  of  the  range  of  an  independent  variable  may 
significantly  bias  the  results  of  a  regression  analysis,  these  points  were  not  taken  into 
account  in  the  analyses  presented  below. 

Two  forms  of  the  regression  equation  relating  the  nondimensional  resistance  in  level  ice, 
R„  to  the  independent  variables,  F  and  C  ,  were  considered.  Form  I  assumed  that  R,  was 
a  second  degree  polynomial  in  Fn 
and  linear  in  Cp.  The  corresponding 
coefficients  are  listed  in  Table  14a 
for  the  total  resistance  and  for  the 
net  resistance.  Because  the  coeffi¬ 
cient  of  F^in  form  I  of  the  regression 
equation  is  relatively  small  and  has 
a  wide  range  of  uncertainty,  as  indi¬ 
cated  by  its  90%  confidence  interval, 
a  second  form  (form  II)  of  the  regres¬ 
sion  equation  was  then  considered 
where  R.is  only  a  linear  function  of 
F2  and  C  ,  i.e.,  R,  =  A  +  B  ■  F2  + 

n  n  n 


Table  13.  Resistance  test  distribution  with  respect 
to  Froude  number  (smooth  model). 


Frouilc  no.  - 

range  CRREL 


Facilili/ 


HSV  A  ISRI  NKK  W  ARC  Total 


0-0.25  6 

0.25-0.50  I 

0.50-0.75  4 

0.75-1.25  5 

1.25- 1.75  1 

1.75-2.25  0 

2.25- 2.75  2 

>  2.75  1 

Tola!  20 


6 

6 

2 

0 

0 

0 

25 


3 

33 


23 

s 

12 

23 

11 

(i 

4 

01 


: 


22 


a.  Data  points  for  full  range  of  Fn. 


Fn 


b.  Form  II  of  regression  equation  for  Fr  <3. 

Figure  13.  Comparison  of  resistance  data  with  smooth  R-class  model  from  all  facilities. 


C  ■  Cn-  The  values  of  the  coefficients  A,  B  and  C  are  listed  in  Table  14b  for  both  the  total  and 
net  resistance. 

From  the  results  of  the  regression  analyses  with  form  I  and  form  II,  it  can  be  seen  that  the 
CRREL  test  results  show  a  somewhat  greater  dependence  of  R.  on  Cn  than  do  those  from 
both  JSRI  and  HSVA.  Consequently,  CRREL  data  yield  a  smaller  value  of  the  intercept  A  of 
both  forms  of  the  regression  equation  as  compared  to  the  intercepts  obtained  for  the  HSVA 
and  JSRI  data. 

The  results  of  these  analyses  are  presented  graphically  for  each  facility  in  Figure  13b, 
where  form  II  of  the  regression  equation  was  plotted  for  the  extreme  values  Cn  =  60  and  Cn 
=  180  of  the  range  of  Cn called  for  in  the  initial  test  program.  In  Figure  14,  the  results  from 
all  the  facilities  are  presented  as  the  measured  values  of  R,  plotted  against  the  values  cal¬ 
culated  by  form  II  of  the  regression  equation  for  the  actual  test  conditions.  The  data  from  all 
facilities  are  in  very  reasonable  agreement  with  one  another,  as  can  be  seen  from  Figure  14 
where  it  is  shown  that  67  data  points  out  of  a  total  of  87  (i.e.,  77%)  fall  within  ±25%  of  the 
predictions. 
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Table  14.  Results  of  regression  analysis  of  resistance  tests  from  all  facilities  (smooth  model;  Fn<2.75). 


a. 

Form  I:  R. 

=  A  +  B  ■ 

F  +  C  F2 

n  n 

+  DCn 

b.  Form  II :  R. : 

=  A  +  B- 

F2  +  C  C 
«  >» 

Laboratory 

A 

B 

c 

D 

r 

Laboratory 

A 

B 

C 

r 

Total  Resistance 

i 

Total  Resistance 

CRREL. 

(90%  C.I.) 

0.066 

±1.425 

0.917 

±1.982 

1.461 

±0.720 

0.0375 

±0.0091 

0.980  i 

! 

CRREL 
(90%  C.I.) 

0.413 

±1.192 

1.777 

±0217 

0.0379 

±0.0089 

0.980 

HSVA 
(90%  C.I.) 

2.070 

±1.090 

-0.002 

±2.929 

2.182 

±1.917 

0.0152 

±0.0071 

0.918 

HSVA 
(90%  C.I.) 

2.070 

±0.712 

2.181 

±0.403 

0.0152 

vO.0069 

0.918 

JSRI 

(90%  C.I.) 

2.685 

±1.367 

0.785 

±1.895 

1.944 

±0.661 

0.0062 

±0.0037 

0.976 

JSRI 

(90%  C.I.) 

3.086 

±0.954 

2209 

±0.163 

0.0060 

±0.0036 

0.977 

NKK 

WARC 

Insufficient  range  for  Cn 

Insufficient  data  range 

i 

All  labs 
(90%  C.I.) 

2.930 

vO.470 

2.122 

±0.121 

0.0089 

±0.0028 

0.959 

All  labs 
(90%  C.I.) 

2.717 

±0.708 

0.450 

±1.115 

1.962 

±0.417 

0.0090 

±0.0028 

0.959 

Net  Resistance 

Net 

Resistance 

CRREL 
(90%  C.I.) 

0.19 

±1.41 

1.09 

±1.96 

0.78 

±0.71 

0.0364 

±0.0090 

0.966 

CRREL 
(90%  C.I.) 

0.60 

±1.19 

1.15 

±0.22 

0.0367 

±0.0089 

0.964 

HSVA 
(90%  C.I.) 

2.04 

±1.06 

0.63 

±2.84 

1.242 

±1.86 

0.0142 

±0.0069 

0582 

HSVA 
(90%  C.I.) 

2.21 

±0.69 

1.64 

±039 

0.0140 

±0.0067 

0.881 

JSRI 

(90%  C.I) 

3.  Go 

±1.11 

0.94 

±154 

1.22 

±054 

0.0042 

±0.0030 

0.968 

JSRI 

(90%  C.I.) 

3.54 

±0.78 

154 

±0.13 

0.0040 

±0.0030 

0.967 

All  labs 
(90%  C.I.) 

2.87 

±0.65 

0.73 

±1.03 

1.23 

±0.417 

0.0073 

±0.0026 

0.971 

All  labs 
(90%  C.I.) 

3.22 

±0.44 

1.49 

±0.11 

0.0072 

±0.0026 

0.933 

Figure  14.  Measured  vs  calculated  total  resistance  (form  II  of  regression  function). 
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Propulsion  tests 

Four  facilities,  namely  CRREL,  HSV  A,  JSRI  and  W  ARC,  carried  out  propulsion  tests  with 
the  smooth  1 :20  scale  R-class  model,  the  results  of  which  were  part  of  the  report  by  the  Ice 
Committee  at  the  17th  ITTC  (ITTC 1984).  Only  the  tests  run  at  HSV  A  were  made  with  a  free- 
running  model;  it  was  connected  to  the  main  carriage  by  only  an  umbilical  cord  for  power 
supply  and  data  transmission.  The  other  three  facilities  carried  out  the  tests  by  the  captive 
model  technique.  At  JSRI  the  propeller  speed  was  adjusted  during  the  tests  to  attain  the  self¬ 
propulsion  point,  i.e.,  so  that  the  net  average  towing  force  or  pull  would  be  zero.  Finally,  the 
tests  at  WARC  were  performed  in  previously  sawn  ice.  The  model  was  towed  by  the 
carriage  at  predicted  speed  and  propeller  rpm  for  the  particular  condi  tions  of  thickness  and 
strength  of  the  parent  ice  sheet,  while  the  pull  exerted  by  the  towing  mechanism,  and 
propeller  thrust  and  torque  were  measured.  The  writers  are  not  familiar  with  the  test 
technique  used  by  WARC  and  are  unsure  whether  the  WARC  data  presented  in  the 
Committee  report  should  be  interpreted  as  being  at  self-propulsion  points  or  not. 

The  data  foi  Iheself-propuision  teats  at  HSVAand  JSRI,  the  published  test  data  of  WARC 
and  the  self-propulsion  test  points  derived  from  the  captive  model  tests  at  CRREL  (as 
explained  in  a  previous  section  of  the  present  report)  are  listed  in  Table  A2.  The  thrust 
deduction  factor,  t,  was  calculated  as 

t  =  l-R/T.  (26) 

P  a 

where  is  the  total  ice  resistance  for  the  test  conditions  of  velocity,  ice  thickness  and  ice 
strength  predicted  using  form  I  of  the  resistance  regression  equation  for  the  resistance  tests 
from  all  facilities.  Ti  is  the  measured  average  thrust. 

Plots  of  the  thrust,  fy  and  torque,  K  ,  coefficients  in  Figures  15a  and  b,  respectively, 
against  the  apparent  advance  coefficient,  /  ,  show  a  linear  relationship: 


•O  • 

— / 


-  o  :j 


JS  p  I 
•  WARC 


a.  Thrust  coefficient,  K;,  vs  apparent  advance  coeffi-  b.  Torque  coefficient,  K  ,  ns  apparent  advance  coefficient, 
dent,}.  J. 


ZCE  “30  r, 

o'  .O  .  _° 

c 

•  ,  o  3  : 


,  R REl  .  jSRI 
m$vA  *  WARC 

_ _ l _ i _ I 

0  ?  0  4 


c.  Thrust  deduction  factor,  t,  t>s  apparent  advance  co¬ 
efficient ,J#. 

Figure  15.  Comparison  of  propulsion  test  results  with  smooth  model, 


(27) 


KfK=A-B  ■/, 


As  seen  from  Figure  15a,  the  thrust  coefficient  measured  in  the  CRREL  tests  appears  to 
be  consistently  greater  than  for  any  other  laboratory,  which  does  confirm  the  suspicion  held 
during  the  tests  that  the  thrust  dynamometers  were  not  functioning  properly.  Similarly, 
from  Figure  15b  it  can  be  seen  that  the  torque  coefficients  measured  at  WARC  are  off  with 
respect  to  the  results  from  the  other  facilities.  Therefore,  the  combined  data  were  analyzed 
with  and  without  those  particular  tests,  and  the  results  listed  in  Table  15  and  shown 
graphically  on  Figure  15. 

The  thrust  deduction  factor,  t,  calculated  according  to  eq  26  is  plotted  against  /ain  Figure 
15c.  No  obvious  variation  of  t  with  /a  is  apparent,  and  for  the  range  of  test  conditions  at  ail 
facilities  the  average  thrust  deduction  factor  is  equal  to  0.168.  Another  approach  is  to  plot 
the  predicted  resistance,  Rp,  against  the  measured  average  total  thrust,  Ta,  as  in  Figure  16, 
which  shows  a  linear  relationship  between  these  two  quantities.  Under  the  assumption  that 


R  +  dR  =  [l-(f  +  dt)]T  (28) 

p  a 

where  dR  is  the  uncertainty  on  Rp  at  the  90%  con¬ 
fidence  interval  and  dt  the  corresponding  uncer¬ 
tainty  on  f,  it  was  found  that  t  =  0.184  ±  0.094. 


T„  (N) 


Figure  16.  Predicted  resistance  of  smooth  model  vs  to¬ 
tal  thrust  measured  during  propulsion  tests. 


Tests  with  roughened  model 

Resistance  tests 

All  six  facilities  (CRREL,  HSVA,  JSRI,  NKK, 
NRCC  and  WARC)  reported  the  results  of  their 
resistance  tests  with  the  roughened  R-class  model 
at  the  1 8th  ITTC,  October  1987,  Kobe,  Japan  (1TTC 
1987).  These  test  results  are  given  in  Table  A3. 
Results  of  regression  analyses  of  the  total  resis¬ 
tance  data  with  forms  I  and  II  of  the  regression 
equation  are  presented  in  Table  16.  The  results  of 
the  analysis  of  the  net  resistance  data  with  form  I 
of  the  regression  equation.  Table  17a,  showed  that 
the  contribution  of  the  term  in  F*  was  negligible, 


Table  15.  Results  of  linear  regres¬ 
sion  analysis  on  K  and  K  (smooth 
R-class  model;  all  laboratories). 


Laboratory 

A 

B 

r 

a.  Thrust  coefficient:  K(  =  A  -  B 

L 

CRREL 
(90%  C.I.) 

0.337 

±0.032 

0.152 

±0.110 

0.624 

HSVA 
(90%  C.I.) 

0.356 

±0.058 

0.408 

±0.159 

0.865 

JSRI 

(90%  C.I.) 

0.324 

±0.018 

0.297 

±0.054 

0.918 

WARC 
(90%  C.I.) 

0.338 

±0.033 

0.338 

±0.096 

0.933 

All  labs 
(90%  C.I.) 

0.341 

±0.020 

0.320 

±0.059 

0.797 

All  labs 
except  CRREL 

0.334 

±0.015 

0.332 

±0.044 

0.904 

b.  Torque  coefficient:  K  ( =  A  -  B 

' 

CRREL 
(90%  C.I.) 

0.0397 

±0.0029 

0.0144 

±0.0102 

0.634 

HSVA 
(90%  C.I.) 

0.0553 

±0.0064 

0.0445 

±0.0175 

0.862 

JSRI 

(90%  C.I.) 

0.0410 

±0.0059 

0.0012 

±0.0174 

0.029 

WARC 
(90%  C.I.) 

0.0537 

±0.0093 

0.0669 

±0.0271 

0.875 

All  labs 
(90%  C.I.) 

0.0448 

±0.0039 

0.0221 

±0.0117 

0.4  lu 

All  labs 
except  WARC 

0.0415 

±0.0034 

0.0080 

±0.0102 

0.206 
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Table  16.  Results  of  regression  analysis  on 
total  resistance  data  with  roughened  model. 


Laboratory 

A 

B 

C 

D 

r 

a.  Form  I:  R,  =  A 

+  B  ■  F 

ft 

+  C  •  F2 

ii 

+  D  C 

CRREL 

2.68 

2.50 

1.00 

0.0174 

0.961 

(90%  C.I.) 

±2.01 

±1.97 

±0.70 

±0.0132 

NRCC 

0.06 

4.93 

0.42 

0.0345 

0.978 

(90%  C.I.) 

±2.53 

±2.88 

±1.03 

±0.0207 

NKK 

3.17 

-0.50 

-2.16 

-0.0247 

0.999 

(90%  C.I.) 

±3.01 

±2.68 

±1.05 

±0.0211 

JSRI 

0.96 

3.44 

0.91 

0.0316 

0.989 

(90%  C.I.) 

±1.27 

±1.41 

±0.51 

±0.0085 

HSVA 

2.19 

1.39 

2.18 

0.0186 

0.987 

(90%  C.I.) 

±2.00 

±2.68 

±0.99 

±0.0128 

WARC 

2.42 

2.54 

3.80 

0.0304 

0.988 

(90%  C.I.) 

±1.69 

±3.37 

±1.98 

±0.0098 

All  except  WARC  1 .41 

2.73 

1.15 

0.0272 

0.966 

(90%  C.I.) 

±1.00 

±1.11 

±0.40 

±0.0064 

b.  Form  11:  R,  =  A  +  B  F^  +  C  C 

*  n  n 


CRREL 
(90%  C.I.) 

4.21 

±1.71 

1.86 

±0.19 

0.0155  — 

±0.0140 

0.954 

NRCC 
(90%  C.I.) 

3.09 

±2.34 

2.12 

±0.35 

0.0300  — 

±0.0176 

0.958 

NKK 
(90%  C.I.) 

2.85 

±2.32 

1.97 

±0.26 

0.0254  — 

±0.0194 

0.984 

JSRI 

(90%  C.I.) 

2.40 

±1.53 

2.11 

±0.19 

0.0339  — 

±0.0116 

0.978 

HSVA 
(90%  C.I.) 

2.84 

±1.54 

2.67 

±0.25 

0.0184  — 

±0.0126 

0.986 

WARC 
(90%  C.I.) 

3.38 

±1.15 

5.25 

±0.42 

0.0277  — 

±0.0093 

0.987 

All  except  WARC 
(90%  C.I.) 

2.87 

±0.87 

2.11 

±0.11 

0.0264  — 

±0.0069 

0.959 

and  a  form  III  of  the  equation,  linear  in  both  Fn  and  Cn,  was  also  applied  and  the  cor¬ 
responding  coefficients  listed  in  Table  17b. 

The  nondimensional  total  resistance  data  are  plotted  in  Figure  17a.  Note  that  the  WARC 
data  points  are  distinctly  higher  than  those  from  all  the  other  laboratories.  In  this  respect, 
the  following  should  be  recognized:  1)  WARC  was  the  first  facility  to  test  the  roughened 
model  and  it  is  possible  that  this  first  series  of  tests  somewhat  smoothed  the  model  so  that 
the  friction  factor  was  greater  at  WARC  than  at  any  other  facility,  and  2)  between  the  tests 
with  the  smooth  model  and  those  with  the  roughened  one,  WARC  changed  its  model  ice 
from  the  columnar  saline  type  to  the  fine-grain  type  and  found  out  that  the  ship  models  had 
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a.  Total  resistance  vs  Froude  number. 


R|t/^Bh^  :  Colculoted 


b.  Measured  vs  calculated  total  resistance. 

Figure  17.  Resistance  tests  with  roughened  R-class  model  at  participating 
facilities. 


to  have  a  lower  friction  coefficient  with  the  new  ice  to  have  the  same  resistance  as  in  the 
previous  model  ice.  For  these  reasons,  W  ARC  data  were  not  included  when  performing  the 
regression  analyses  on  all  available  test  data.  The  corresponding  regression  line  (form  I  of 
the  equation)  is  also  shown  in  Figure  17a.  Figure  17b  presents  the  comparison  between  the 
resistance  given  by  the  regression  equation  derived  and  that  measured  at  all  facilities  except 
WARC.  It  can  be  seen  that  all  data  points  fall  within  ±  25%  of  the  calculated  values. 

Propulsion  tests 

Only  three  facilities,  namely  CRREL,  JSRI  and  NRCC,  have  reported  results  of  propul¬ 
sion  tests  with  the  roughened  model,  the  results  of  which  are  listed  in  Table  A4. 

The  thrust  and  torque  coefficients,  K(  and  K  ,  are  plotted  against  the  apparent  advance 
coefficient  Ja  on  Figure  18,  and  the  results  of  tire  linear  regression  analyses  of  and 
according  to  eq  28  are  listed  in  Table  18.  The  regression  line  for  all  laboratories  is  shown  on 
Figure  18.  It  can  be  observed  that  the  results  on  and  K  obtained  at  the  various  facilities 
are  in  excellent  agreement  with  one  another,  and  companson  with  the  data  obtained  with 
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Table  17.  Results  of  regression  analysis  on  net 
resistance  data  with  roughened  model. 


Laboratory 

A 

B 

C 

D 

r 

a.  Form  I :  R 

:.  =  A  +  B  •  Fn 

+  C  F„2 

♦  D  C, 

CRREL 

2.74 

2.91 

072 

0.0157 

0.944 

(90%  Cl.) 

±1.64 

±1.61 

±0.57  ±0.0108 

NRCC 

0.32 

579 

-0.30 

0.0314 

0.966 

(90%  C.I.) 

±2.37 

±2.70 

±0.097  ±0.0130 

NKK 

4.06 

-0.15 

152 

0.0162 

0.976 

(90%  Cl.) 

±2.75 

±2.45 

±0.96  ±0.0192 

JSRI 

1.26 

4.09 

0.06 

0.0272 

0.990 

(90%  C.I.) 

±0.92 

±1.02 

±0.37  ±0.0062 

HSVA 

2.32 

1.75 

1.49 

0.0162 

0.987 

(90%  Cl.) 

±1.62 

±2.18 

±0.80  ±0.0104 

WARC 

2.51 

3.24 

2.80 

0.0285 

0.988 

(90%  Cl.) 

±1.49 

±2.97 

±1.74  ±0.0086 

All  except  WARC 

1.59 

3.19 

0.38 

0.0244 

0.946 

(90%  C.I.) 

±0.93 

±1.03 

±0.38  ±0.0059 

b.  Form  III:  R. 

=  A  +  B 

'  F„  +  C 

c 

n 

CRREL 

2.41 

3.50 

0.0162 

_ 

0.943 

(90%  C.I.) 

±1.37 

±0.41 

±0.0106 

NRCC 

0.77 

4.49 

0.0306 

_ 

0.965 

(90%  C.I.) 

±1.79 

±0.68 

±0.0122 

NKK 

1.48 

3.57 

0.0244 

_ 

0.949 

(90%  C.I.) 

±2.99 

±0.89 

±0.0249 

JSRI 

170 

4.25 

0.0271 

— 

0.990 

(90%  Cl.) 

±0.82 

±076 

±0.0060 

HSVA 

051 

5.66 

0.0182 

_ 

0.973 

(90%  Cl.) 

±1.76 

±0.76 

±0.0140 

WARC 

0.83 

7.91 

0.0337 

_ 

0.981 

(90%  C.I.) 

±1.29 

±0.75 

±0.0097 

All  except  WARC 

1.08 

4.21 

0.0249 

— 

0.944 

(90%  C.I.) 

±0.79 

±0.27 

±0.0060 

the  smooth  model  shows  no  significant  effect  of  the  hull  roughness  on  the  thrust  and  torque 
coefficients. 

The  thrust  deduction  factor,  t,  was  calculated  by  eq  26  where  Rpis  the  total  ice  resistance 
given  by  form  I  of  the  resistance  equation,  namely 

R  /yBtf=  1.413  +  2.73  F  +  1.149  F2  +  0.0272  C  (29) 

p  '  i  n  n  n 

and  is  plotted  versus  /  on  Figure  19.  Considering  the  inherent  uncertainty  in  values  of  t  for 
tests  in  ice,  there  is  remarkable  agreement  between  the  three  laboratories.  In  contrast  with 
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a.  Thrust  coefficient,  K(/  vs  apparent  advance  coeffi¬ 
cient  ,]a. 


b.  Torque  coefficient,  K^,  vs  apparent  advance  coeffi¬ 
cient,  Ja. 


Figure  18.  Comparison  of  propeller  coefficients  from  propulsion  test  results  with  roughened  model. 


the  results  for  the  smooth  model,  the  scatter  in  t  is  relatively  moderate  and,  in  particular, 
there  are  no  negative  values.  There  is  no  evident  effect  of  /  on  t,  which  can  be  considered 
constant,  equal  to  0.182  with  a  standard  deviation  of  0.073.  This  is  identical  to  the  result  of 
t  =  0.184  ±  0.094  obtained  for  the  smooth  model. 

Another  way  of  determining  t  is  to  plot  the  total  resistance  predicted  by  eq  29  against  the 
corresponding  measured  total  thrust,  T  ,  as  in  Figure  19b.  Linear  regression  of  these  data 
yielded  f  =  0.178  with  a  90%  confidence  interval  of  ±0.021 ,  equivalent  to  the  previous  result. 


Comparison  with  full-scale  data 

As  noted  earlier,  the  values  of  the  propeller  co¬ 
efficients,  /C  and  K  ,  obtained  from  model  tests  with 
both  the  smooth  and  roughened  models  were  nearly 
identical.  In  Figure  20,  these  model  test  values  are 
compared  with  full-scale  data  calculated  from  the 
measurements  presented  by  Michailidis  and  Murdey 
(1981)  and  with  results  from  model  tests  in  clear  water 
(Murdey  1980).  It  can  be  seen  that  the  values  of  Kt 
obtained  in  both  types  of  model  tests,  i.e.,  in  level  ice 
and  in  clear  water,  are  practically  identical  and  match 
well  the  full-scale  values,  especially  for  /  >  0.25,  ap¬ 
proximately.  For  /  <  0.2,  the  full  scale  values  of  K(are 
somewhat  lower  than  those  obtained  in  model  tests. 
On  the  other  hand,  the  values  of  the  torque  coefficient, 
K ' ,  measured  in  model  tests  in  level  ice  are  consis¬ 
tently  greater  than  those  measured  in  clear  water 
tests,  and  the  latter  are  in  excellent  agreement  with  the 

full  scale  values  of  K  . 

q 

T  he  ship  performance  at  full  power  predicted  from 
model  test  results  are  presented  in  Figure  21  as  maxi¬ 
mum  velocity  versus  ice  thickness  for  an  assumed 
value  of  450  kPa  of  the  ice  strength,  the  average  of  a. 
reported  by  Michailidis  and  Murdey.  Figure  21  also 
shows  the  full-scale  data  reported  by  Michailidis  and 
Murdey  and  those  presented  in  the  1987  ITTC  Ice 
Committee  report  (ITTC  1987)  based  on  all  available 


Table  18.  Results  of  linear  re- 

f;ression  analysis  on  and  K 
roughened  R-class  model;  alt 
facilities). 


Laboratory 

A 

B 

r 

a.  Thrust  coefficient:  K(  =  A  -  B 

I 

CRREL 
(90%  C.I.) 

0.358 

±0.037 

0.365 

±0.117 

0.965 

JSRI 

(90%  C.I.) 

0.354 

±0.016 

0.326 

±0.051 

0.974 

NRCC 
(90%  C.I.) 

0.304 

±0.038 

0.232 

±0.109 

0.929 

Aii  labs 
(90%  C.I.) 

0.344 

±0.015 

0.319 

±0.047 

0.936 

b.  Torque  coefficient:  =  A  -  B 

Ja 

CRREL 
(90%  C.I.) 

0.0402 

±0.0052 

0.0026 

±0.0166 

0.178 

JSRI 

(90%  C.I.) 

0.0425 

±0.0075 

0.0031 

±0.0234 

0.087 

NRCC 
(90%  C.I.) 

0.0484 

±0.0102 

0.0211 

±0.0295 

0.644 

All  labs 
(90%  C.I.) 

0.0432 

±0.0040 

0.0072 

±0.0124 

0.220 

30 


a.  Calculated  In/  cej  26.  />.  From  plot  of  R^versus  T(. 

Figure  19.  Thrust  deduction  factor. 


a.  Thrust  coefficient,  Kf,  t»s  apparent  advance  coeffi-  b.  Torque  coefficient,  K  ,  us  apparent  advance  coefficient, 
dent,).  T,  q 

Figure  20.  Comparison  uetween  model  test  results  and  full-scale  data. 


h,  (m ) 


a.  Based  on  smooth  model  test  results.  b.  Based  on  roughened  model  test  results. 

Figure  21.  Comparison  between  predicted  performance  in  leivl  ice  and  full-scale  measurements  (a  = 
450  kPa). 


reports  of  field  trials  with  the  R-class  icebreakers.  It  should  be  noted  that  the  thickness 
interval  indicated  by  a  solid  line  segment  on  Figure  21  for  the  data  of  Michailidis  and 
Murdey  represents  the  thickness  of  the  snow  layer  overlaying  the  sheet  ice.  The  ship 
performance  predictions  were  calculated  from  the  model  test  results  with  both  the  smooth 
(Fig.  21a)  and  rough  (Fig.  21b)  models,  and  using  the  values  of  K ^  and  t  obtained  from 
model  test  in  both  level  ice  and  clear  water.  From  Figure  21,  it  can  be  concluded  that: 

1.  The  predictions  based  on  the  resistance  tests  results  with  the  rough  model  match  the 
full-scale  trial  data  better  than  those  based  on  resistance  tests  results  with  the  smooth  model. 

2.  For  ice  thickness  less  than  0.7  m,  the  "rough  model"  predictions  based  on  propeller 
coefficients  from  tests  in  clear  water  are  in  near-perfect  agreement  with  the  full-scale 
measurements.  As  the  ice  thickness  increases  beyond  0.7  m,  and  the  maximum  ship  speed 
at  full  power  decreases  accordingly,  the  full-scale  data  agree  more  and  more  with  the  pre¬ 
dictions  based  on  propeller  coefficients  from  model  tests  in  level  ice.  This  observation 
confirms  the  concern  that  model  propulsion  tests  in  level  ice  performed  in  relatively  thin 
level  ice,  or  at  high  speed,  yield  propeller  coefficients  that  are  too  pessimistic. 

3.  For  the  field  conditions  documented  on  Figure  21,  the  effect  of  snow  on  the  ship 
performance  can  be  accounted  for  by  increasing  the  ice  thickness  by  half  the  snow  thickness. 


CONCLUSIONS  AND  RECOMMENDATIONS 
CRREL  tests 

Following  resistance  and  propulsion  tests  in  level  ice  with  a  smooth  (friction  factor  of 
about0.04)  model  of  the  R-class  icebreaker,  similar  tests  were  performed  with  a  roughened 
model.  From  the  results  of  these  tests  the  following  observations  can  be  made. 

1.  The  ice  friction  tests  indicated  that  the  apparent  friction  coefficient,/^  defined  as  the 
ratio  of  the  tangential  force,  T,  to  the  normal  load,  N,  decreases  with  increasing  N.  In  other 
words,  the  tangential  force  is  not  directly  proportional  to  the  normal  load  but  a  linear 
function  of  it.  The  intercept,  Tq,  was  interpreted  as  an  adher  n  stress,  the  exact  origin  of 
which  remains  to  be  clarified,  while  the  slope  was  interpreted  as  the  actual  friction 
coefficient.  The  friction  factor  of  the  original  "smooth"  model  had  been  measured  at  0.04; 
that  of  the  roughened  model  was  approximately  0.07. 

2.  The  increase  in  nondimensional  resistance  due  to  the  increase  in  the  friction  coefficient 
decreased  with  increasing  nondimensional  ice  strength,  Cn.  For  Cn  =  60,  the  resistance  of 
the  roughened  model  was  10  to  20%  greater  than  that  of  the  smooth  model,  while  for  Cn  = 
180  the  resistance  of  both  models  was  nearly  the  same. 

3.  The  propeller  characteristics,  fy  Kqand  f,  of  both  models  were  practically  identical, 
showing  no  effect  of  the  increase  in  hull  roughness. 

Comparative  studies 

Both  models  had  been  tested  by  several  research  facilities  both  in  resistance  and  in 
propulsion.  Comparison  of  the  test  results  between  these  facilities  can  be  summarized  as 
follows. 

1 .  For  the  most  part,  the  resistance  and  propulsion  test  results  from  any  facility  fall  within 
±25%  of  the  overall  average  test  results.  Furthermore,  for  given  test  conditions  of  ice 
thickness,  ice  strength  and  ship  speed,  the  variation  in  test  results  between  research  facilities 
is  of  the  same  order  as  that  at  any  one  facility. 

2.  The  model  propulsion  test  results  of  all  facilities  showed  that  the  propeller  thrust 
coefficient  in  level  ice  was  practically  the  same  as  in  clear  water,  while  the  torque  coefficient 
in  level  ice  was  nearly  constant  and  significantly  greater  than  in  clear  water.  Significant 
uncertainty  in  the  thrust  deduction  factor,  especially  witi  i  thesmooth  model,  deduced  from 
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i 

j  the  results  of  propulsion  tests  in  level  ice  has  to  be  attributed  primarily  to  unavoidable 

|  uncertainty  in  the  model  resistance  test  results.  Overall,  the  thrust  deduction  factor  in  level 

\  ice  was  found  to  be  approximately  constant,  equal  to  0.1 8  ±0.07,  rather  than  increasing  with 

ship  speed  (or  apparent  advance  coefficient)  as  in  clear  water. 

!  3.  Ship  performance  prediction  presented  as  maximum  ship  speed  versus  ice  thickness 

based  on  roughened  model  resistance  and  propeller  characteristics  in  clear  water  showed 
excellent  agreement  with  available  full-scale  trials  data  for  ice  thickness  up  to  0.7  m.  For 
■  thicker  ice,  good  comparisons  were  obtained  when  the  propeller  characteristics  from  tests 

in  ice  were  used  in  the  prediction  algorithms.  These  observations  would  seem  to  indicate 
that,  compared  to  full  scale,  ice  effect  on  propulsion  is  exaggerated  in  model  tests  at  rel¬ 
atively  high  speed  in  ice  of  low  to  medium  thickness,  but  is  correctly  modeled  in  tests  in 
thick  ice. 

It  is  extremely  encouraging  to  find  that  ship  testing  in  ice  has  reached  a  level  of  reliability 
where  different  facilities,  using  somewhat  different  model  ice  and  slightly  different  testing 
techniques,  obtain  overall  results  that  are  comparable,  both  in  resistance  and  propulsion, 
with  an  uncertainty  in  the  results  that,  in  the  opinion  of  the  writers,  is  inherent  to  the 
icebreaking  phenomenon.  However,  progress  in  the  development  of  model  ice,  modeling 
techniques,  and  data  analysis  and  interpretation  is  still  needed,  especially  concerning 
ice-propeller  interaction,  which  is  usually  exaggerated  at  the  model  scale.  This  area  of 
research  is  one  that  is  and  ought  to  be  receiving  further  attention  from  all  ice  research 
facilities  in  order  to  achieve  a  level  of  confidence  in  the  results  of  model  tests  in  ice  and  their 
extrapolation  to  full  scale  comparable  to  that  enjoyed  by  model  tests  in  a  traditional,  open- 
water  towing  tank. 
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APPENDIX  A:  TEST  RESULTS  (ALL  LABORATORIES) 


Table  A2.  Self-propulsion  points  of  smooth  model, 
a.  Dimensional  data. 
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Table  A2  (cont'd).  Self-propulsion  points  of  smooth  model 
b.  Nondimensional  data. 


Table  A3  (cont'd).  Resistance  test  results  with  roughnened  model 
b.  Nondimensional  data. 
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Table  A4.  Self-propulsion  points  of  roughnened  model, 
a.  Dimensional  data. 
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b.  Nondimensional  data. 
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